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Abstract 


Growing energy crops on marginal land has been promoted as a way of ensuring that biomass production 
involves an acceptable and sustainable use of land. Saline and saline-prone agricultural lands represent an 
opportunity for growing energy crops avoiding the displacement of food production and contributing to restora- 
tion of degraded land. Giant reed (Arundo donax L.) is a perennial grass that has been proposed as a promising 
energy crop for lignocellulosic biomass production while its tolerance to salinity has been proved. In this work, 
the identification of surplus saline lands that could be irrigated with saline waters for growing tolerant-energy 
crops (giant reed) in the mainland of Spain and the assessment of the agronomically attainable yield in these 
limiting growing conditions were undertaken. To this purpose, a GIS analysis was conducted using geo- 
databases related to saline areas, agro-climatic conditions, irrigation water requirements, agricultural land 
availability, restrictions regarding the range of electrical conductivity tolerated by the crop, competition with 
agro-food crops and irrigation water provisions. According to the approach developed, the irrigated and saline 
agricultural area available and suitable for biomass production from giant reed amounted up to 34 412 ha. The 
agronomically attainable yield in these limiting conditions was estimated at 12.7 — 22.2 t dm ha’! yr and the 
potential production of lignocellulosic biomass, 597 338 t dm yr’. The methodology followed in this study can 
be applied to other target regions; it allows the identification of this type of marginal lands, where salinity-toler- 
ant plant species could be grown for bioenergy purposes, avoiding competition with agro-food crops, and where 
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soil restoration measurements should be undertaken. 
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Introduction 


Growing energy crops on marginal land has been pro- 
moted as a way of ensuring that biomass production 
involves an acceptable and sustainable use of land. The 
condition of ‘marginal’ may be ambiguous and site- 
specific; a piece of land may be marginal in a specific 
region or for a certain purpose, but it may be suitable 
for a different objective and even its quality can be con- 
sidered no marginal in another location. Therefore, 
availability and suitability of marginal land for biomass 
production is surrounded by uncertainties among the 
wide-ranging estimates (Dale et al., 2010; Lewis & Kelly, 
2014). Actually, there are three separate definitions of 
marginal land: land unsuitable for food production 
(‘because the land is not productive enough’); ambigu- 
ous low-quality land (‘land not necessarily unsuitable 
for food production but where food production is less 
productive’); and economically marginal land (‘land 
where cost-effective agricultural production is not 
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possible under a given set of conditions’) (Turley et al., 
2010; Shortall, 2013). Even though there is not a 
complete agreement on the definition of marginality, 
saline and saline-prone agricultural lands are certainly 
covered by the former concepts and they may represent 
an opportunity for growing energy crops avoiding the 
displacement of food production and contributing to 
restoration of degraded land. 

The causes of soil salinization may be divided into 
natural and human-induced factors. Natural causes are 
subdivided into land (climate, geology, hydrology, etc.) 
and soil factors (texture, structure, compaction rate, etc.) 
(Toth et al., 2008). The latter factors are mainly based on 
land management and degradation. Thus, high salt con- 
tents (sodium, potassium, magnesium and calcium, 
chloride, sulphate, carbonate) in parent materials or 
groundwater lead to primary salinization (i.e. by natural 
processes) whereas human practices such as inappropri- 
ate irrigation or overexploitation of aquifers drive to 
secondary salinization (Toth et al., 2008). 

Some ways for remediating saline soils and for tak- 
ing advantage of them are proper agronomic practices, 
use of salt-tolerant crops and phytoremediation 
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(Hasanuzzaman et al., 2014). Nevertheless, the response 
of plants to water and saline stress conditions may be 
the closure of stomata to reduce water losses, which 
would inhibit photosynthesis. This effect leads to reduc- 
tion in plant growth and productivity. In addition, salt 
accumulation in soils reduces its hydraulic conductivity, 
which diminishes soil water availability (Cosentino 
et al., 2013). The development of biomass production 
from crops tolerant to limiting conditions — particularly 
drought and salinity — is among the main targets for 
crop management in Mediterranean areas (Araus et al., 
2002; Cosentino et al., 2013). 

Giant reed (Arundo donax L.) is a perennial grass that 
has been proposed as a promising energy crop for ligno- 
cellulosic biomass production in marginal, abandoned or 
set-aside agricultural land while its tolerance to salinity 
has been evidenced within the EU FP7 OPTIMA project 
(Curt et al., 2014; Sanchez et al., 2015). Actually, the toler- 
ance of giant reed to salinity conditions and its physio- 
logical effects have been tested in several studies and 
regions, such as Italy, Spain, USA and Australia (Wil- 
liams et al., 2009; Zema et al., 2012; Cosentino et al., 2013; 
Sanchez et al., 2015). In the study by Williams et al. 
(2009), it was found that some robust giant reed stands 
were able to tolerate electrical conductivities of soil 
extracts of 35-45 dS m™!. Other studies (Calheiros et al., 
2012) have showed the potential of giant reed for polish- 
ing organic matter and nitrogen from a_high-salinity 
effluent from a tannery’s wastewater (13.3-19.3 dS m~’) 
in constructed wetlands. Moreover, giant reed has been 
recently cited in the European Union Directive 2015/ 
1513 as ‘non-food cellulosic material’; therefore, its 
cultivation for biofuels can be counted for the national 
targets of the EU Member States and its contribution is 
not subjected to the 7% limit of the final consumption of 
energy in transport in 2020. 

Regarding the study of bioenergy potential produc- 
tion in degraded lands, different methodologies and 
assessments at different World regions and scales may 
be found in the bibliography (Milbrandt & Overend, 
2009; Wicke et al., 2011; Zhuang et al., 2011; Nijsen et al., 
2012; Tenerelli & Carver, 2012; Kang etal., 2013). 
Methodologies for identification of saline soils at regio- 
nal scale are also quite abundant in the literature, 
mainly based on geographic information systems and 
remote sensing tools (Douaoui et al., 2006; Acosta et al., 
2011; Bouaziz et al., 2011; Melendez-Pastor et al., 2012; 
Ochieng ef al., 2013). On the other hand, the specific 
case of soil salinity assessment for biomass production 
from dedicated energy crops is less widespread: for 
example, a worldwide analysis of the extent and loca- 
tion of salt-affected soils and its potential to produce 
biomass from forestry plantations was conducted by 
Wicke et al., 2011. To the best of our knowledge, no 


assessment of surplus saline soils for biomass produc- 
tion has been conducted at a national scale in Spain so 
far. By surplus saline lands, it is meant those areas 
classified simultaneously as saline and saline-prone 
lands that would not represent competition with the 
agro-food production if giant reed were grown. 

The specific objectives of this work were the follow- 
ing: (i) to develop a method to identify and assess the 
saline or saline-prone marginal lands that could be irri- 
gated with saline waters for growing tolerant-energy 
crops such as giant reed, and (ii) to apply this method 
to the case study of Spain in order to estimate the ligno- 
cellulosic biomass potential production of giant reed. 


Materials and methods 


Identification of saline and saline-prone areas 


Saline agricultural lands (saline soils) and saline-prone lands 
(agricultural land irrigated with brackish or saline waters) were 
assessed altogether in order to estimate the area of land sub- 
jected to saline degradation processes that might be suitable for 
growing giant reed. 

Regarding saline soils, the database of saline and sodic soils 
in the European Union (Toth et al., 2008) was compiled and 
polygons referred to mainland Spain were extracted (Fig. 1). 
This database displays the calculated percentage share of salt- 
affected areas of the polygons within the Soil Mapping Units 
from the European Soil Database (European Commission and 
the European Soil Bureau Network, 2004). The three classes of 
soil salinity (more than 50% of the area with saline soils, less 
than 50% and potentially salt-affected areas) were selected as 
potential saline lands. 

The identification of saline-prone land was based on the 
measurements taken by the control networks of water quality 
established by each of the 15 River Basin Districts (RBD) in 
Spain (Fig. 2) for the accomplishment of the Water Framework 
Directive for water protection and management (Directive 
2000/60/CE). 

The above-mentioned networks are divided in two sections, 
superficial and groundwater water control networks. In both 
cases, they evaluate different parameters of water quality, 
including water electrical conductivity (ECw) measured at 
25 °C. Results of these measurements from some RBD control 
networks are available on their respective websites, but some 
databases need to be officially requested. Thus, data of ECw, 
measured at 25 °C, were gathered and compiled for the time 
period 2010-2013 from each RBD for both surface and ground- 
water control networks. The annual mean values for 2010 were 
calculated, geo-referenced and depicted in a geographic infor- 
mation system (GIS) environment (ArcGIS — Environmental 
System Research Institute — ESRI™). Average values of 2010 
were taken for this work as the 2010 data set was the largest 
and most complete one from the compiled database for the 
analysed time period. These values are shown in Fig. 3. 

Geostatistical analyses were conducted in order to interpolate 
ECw point data to continuous data in raster layers using the 
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Fig. 1 Map of saline soils in mainland Spain. Geo-referenced information extracted from the European Soil Database (Toth et al., 


2008) 
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Fig. 2 River Basin Districts in mainland Spain. 


Geostatistical wizard (ESRI™). Databases for surface and ground- 
water control networks were interpolated separately but also 
jointly with the integrated ECw database. ECw data from 3618 
and 2242 points were collected from the groundwater and sur- 
face control networks, respectively. In both cases, 1% of the total 
point data were randomly selected and laid away as a subset for 
the subsequent model validation. The interpolation methods 


tested were deterministic (Inverse Distance Weighted, Spline, 
Local Polynomial and Radial Basis Function) and geostatistical 
(Krigging). Mean relative errors (MRE) and root-mean-squared 
errors (RMSE) of each interpolation method were calculated. 
Once an ECw continuous data set was obtained, the average 
value in each area of irrigated land was calculated by means of 
the ‘zonal statistic’ tool, with the aim to assign values for ECw 
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Fig. 3 Mean electrical conductivity of water (ECw) according to the control networks of water quality. 


to each polygon. For this purpose, the class ‘permanently 
irrigated agricultural land’ from the Corine Land Cover 2006 
vector data (http://www.eea.europa.eu/data-and-maps/data/ 
cle-2006-vector-data-version-3) was used (henceforth, CLC2006). 
A correlation between the irrigation water origin: surface irri- 
gation, groundwater irrigation and undifferentiated irrigation 
(a combination of the two types) and water quality control 
network was assumed. Thus, according to the information 
gathered from the Agricultural Demand Units (ADU) — see sec- 
tion Analysis of resource availability — one type of water origin 
was assigned to each irrigated land polygon and therefore its 
corresponding ECw. 

Then, the electrical conductivity of saturated soil extract 
(ECe) was calculated from ECw in each polygon of irrigated 
land, by means of the equation by Ayers & Westcot (1985). 


ECe = ECw x 1.5 (1) 


According to consulted bibliography, several clones of giant 
reed are tolerant to salinity and are able to grow with ECe val- 
ues up to 8-9 dS m™! (Cosentino et al., 2013), whereas irrigated 
cereals — wheat, barley, triticale — can tolerate 6 dS m™! at the 
most (Rhoades et al., 1992). 

Therefore, the saline agricultural land available for growing 
giant reed was estimated by the selection of those polygons 
with an estimated ECe between 6 and 9 dS m | which cen- 
troids (geometric centres) were within the polygons of saline 
soils mapped by Toth et al. (2008). 


Estimation of crop requirements 


Giant reed is known to be a perennial species adapted to limit- 
ing growing conditions, but it also needs some basic resources 


to be able to thrive and produce enough biomass. The method- 
ology developed for this section aimed at the analysis and esti- 
mation of the crop requirements, in terms of crop cycle length 
and irrigation needs of giant reed in mainland Spain, to 
identify the irrigated land suitable for its cultivation. 

The estimation of the crop cycle length was based on the cal- 
culation of the growing degree days (GDD, °C day), commonly 
used for measuring crop phenology and development with the 
concept of heat units (McMaster & Wilhelm, 1997) needed for 
plant growth. The irrigation needs were estimated from the 
crop coefficient (Kc), defined as the ratio of evapotranspiration 
(ETc) from a specific crop or soil surface to a reference ET 
value (Allen et al., 2005). 

The input data for the estimation of the crop cycle length 
and Kc of giant reed in Spain were extracted from experiments 
conducted in the Mediterranean region. Kc and water-use effi- 
ciency of giant reed for different growth stages were estimated 
from data of a field experiment conducted in Pisa in 2010 and 
2011 (Triana et al., 2014). The growth stages were identified as: 
initial (from crop sprouting to the beginning of stem elonga- 
tion), crop development (stem elongation), mid-season (from 
the end of stem elongation to the beginning of canopy senes- 
cence) and late season (from canopy senescence to the end of 
water uptake). In our work, GDDs for each stage were calcu- 
lated from data of Pisa-S.Giusto weather station (www. 
tutiempo.net/clima) assuming 10 °C as base temperature (Di 
Nasso et al., 2011). Average data of both years were calculated 
and extrapolated to Spanish conditions (Table 1). 

Following this approach, the length of giant reed growth 
cycle in Spain was spatially determined either by GGD or stage 
duration (days). For GGD calculation, the mean temperature 
for each set of 10 days was interpolated from mean monthly 
temperatures raster data set (National Meteorological Agency 
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Table 1 Growth stages and crop coefficient (Kc) (+ standard deviation) of giant reed reported by Triana et al. (2014) and Growing 


Degree Days (GDD) calculated from climate data for this work 


Year Growth stage Starting date 
2010 Initial 29-Mar 
Crop development 12-May 
Mid-season 06-Jun 
Late season 17-Sep 
Total 21-Oct 
2011 Initial 20-Mar 
Crop development 21-May 
Mid-season 12-Jun 
Late season 06-Sep 
Total 21-Oct 
Average Initial 24-Mar 
Crop development 17-May 
Mid- season 09-Jun 
Late season 12-Sep 
Total 21-Oct 


of Spain — AEMET, 2010a) using Gommes’ interpolation 
method (2) (Gommes, 1983): 


Ti0-1 = (5 x TmM-1 +26 x Tm-2 —4x Tm-3)/27 
T10-2 = (—Tm-1 +29 x Tm-2 = Tm_3)/27 (2) 
T10-3 = (—4 x Tm-1 +26 x Tm-2 +5 x Tm-3)/27 


where Tio-1; Ti9-2 and Tio-3 are the temperatures in each set of 
10 days within a month and Tyy-1, Tm-2 and Ty-3 are the mean 
monthly temperatures in the corresponding 3-month interval. 
The starting date of the crop cycle was considered to be 
dependent on the mean minimum temperature (National Mete- 
orological Agency of Spain — AEMET, 2010b), and the limit 
was set at 3 °C following the criterion by Emberger (1955) 
when frost period is less frequent. Subsequently, GDD and the 
Kc for each set of 10 days were geo-spatially estimated. 
Irrigation needs for agriculture depend on both crop require- 
ments and climatic characteristics of the studied area. Hence, 
the monthly effective precipitation — Pef — (fraction of rainfall 
that is really available for plants) was estimated following 
Brouwer and Heibloem method (3) (Dastane, 1978) from 
monthly accumulated data (P) provided in a raster data set 
(National Meteorological Agency of Spain - AEMET, 2010c): 


Pef =0.8x P—25 if: P>75mmmonth ! 


(3) 
Pef = 0.6 x P—10 


if P <75mmmonth | 

In addition, monthly evapotranspiration (PET) was also com- 
piled from the data sets available at the Consortium for Spatial 
Information (CGIAR-CSI, http://www.cgiar-csi.org/) which 
are based on a globally modelled evapotranspiration (Trabucco 
et al., 2008) from a spatialized implementation of the Harg- 
reaves (Hargreaves, 1994) evapotranspiration equation (Zomer 
et al., 2008). Only the monthly PET values for the giant reed 
crop cycle were gathered. 

Gommes’ interpolation method for accumulating variables 
was then used for calculating Pef and PET values for each set 
of 10 days from monthly values (Gommes, 1983): 


Ke Stage duration (days) GDD (°C day) 
0.67 + 0.07 44 154 
0.91 + 0.14 25 182 
1.93 + 0.09 103 1320 
1.38 + 0.13 34 244 
- 206 1900 
0.41 0.06 43 262 
1.11 + 0.14 22 239 
1.55 + 0.14 86 1158 
0.98 + 0.12 45 427 
- 196 2086 
0.54 + 0.18 44 208 
1.01 + 0.14 24 211 
1.74 + 0.27 95 1239 
1.18 + 0.28 40 336 
- 201 1993 


PET}0-1 = (5 x PETy-1 +26 x PETy_2 —4x PETm-3)/81 
PET 0-2 = (—PETq_1 + 29 x PETm-2 — PETw_3) /81 
PET j0_3 = (-4 x PETy-1 +26 x PETy-2 +5x PETw-3)/81 


(4) 


where PETjo-1, PETi0-2 and PET10-3 represent the evapotranspi- 
ration values in each set of 10 days within a month and PETy- 
1, PETm-2 and PETy,3 are the mean monthly evapotranspiration 
in the corresponding 3-month interval. 

Then, the irrigation needs of giant reed per set of 10 days 
were calculated according to the expressions (5) and (6). 


ETc = ETo x Ke (5) 


Ar = ETc — Pef (6) 


where ETc is the crop evapotranspiration and Ar represents 
crop irrigation requirements for an optimal yield. 

Water irrigation requirements for each set of 10 days were 
then calculated using ESRI Spatial Analyst tool. Finally, the 
zonal statistics tool was used for the calculation of the mean 
values of irrigation requirements in each polygon of irrigated 
agricultural land during the whole crop cycle. 


Analysis of resource availability 


This section aims at the compilation of geo-referenced data 
regarding the basic agrarian resources for the cultivation of 
giant reed in Spain: land and water. 

The irrigated agricultural land was based in the Corine Land 
Cover 2006 vector data. The class ‘permanently irrigated agri- 
cultural land’ was considered as the agricultural land poten- 
tially available for crops with accessible irrigation systems. 

As regard to availability of water for crop irrigation, each 
RBD in Spain is in charge of elaborating its respective Hydro- 
logical Plan (2010-2015) (HP) for water management, in accor- 
dance with the Water Framework Directive of the European 
Union. Within each Plan, water demand is estimated for 
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human consumption, agricultural (crops + livestock) uses, 
industrial uses and ‘others uses’. According to the water 
demand for the different uses and the water availability in the 
river basin, each RBD designates the amount of water for each 
of the above-mentioned uses, which includes the water provi- 
sion for crops irrigation. 

The information provided by the Hydrological Plans in 
Spain is rather heterogeneous: some Plans specifically allocate 
the water provision for each of the existing crops in the river 
basin (considering the water requirements of each crop); others 
specify the water provision per type of irrigation system; 
certain Plans differentiate between surface and groundwater 
provision; water provisions refer either gross and net amount 
depending on the analysed Plan; the geo-referenced unit varies 
and the water provision is allocated per Agrarian County (an 
ensemble of municipalities) or per Agricultural Demand Unit — 
ADU (irrigation area with similar characteristics in terms of 
location, water origin or irrigation community to which the 
area belongs). Hence, a homogenization process was conducted 
in this work to build a geo-referenced database of general 
water provision for crop irrigation at national scale. Assump- 
tions made for this purpose were the following: 


e When the HP made the distinction between irrigation 
systems, sprinkling-type irrigation was chosen. 


e Linked to the type of irrigation system, 75% efficiency was 
assumed (Confederación Hidrográfica del Duero, 2012) for 
sprinkling-type irrigation in case the water provision was 
given in gross figures. 


e When no differentiation between the origins of the 
irrigation water (surface or groundwater), undifferentiated 


irrigation was assigned and consequently, the ECw corre- 
sponding to the integrated database as well — see section 
Identification of saline and saline-prone areas. 


This way the water provision for crops irrigation in each 
Hydrological Plan in Spain was gathered, homogenized and 
introduced in the GIS environment (Fig. 4). 

As the geo-referenced units — to which water provisions 
were allocated — slightly differed from the polygons of irrigated 
agricultural land from CLC2006 data set, an assignation pro- 
cess between the two layers was conducted. To this purpose, 
water provision for the three types of water origin (surface, 
groundwater and undifferentiated) was converted into continu- 
ous raster data sets with a pixel size of 100 m. The zonal statis- 
tics tool was used to calculate the mean water provision for 
each polygon of irrigated land in CLC2006 as well as the 
number of pixels from each type of water origin within each 
polygon. In those with more than one water origin type, the 
predominant one was assigned. 


Estimation of biomass potential production 


Once the crop requirements were assessed in terms of crop 
cycle length and irrigation water and the availability of the 
related agricultural resources was evaluated, the potential 
production of biomass (stems and leaves) under limiting condi- 
tions was estimated. The methodology consisted in the applica- 
tion of yield constraints on the crop potential yield, similarly to 
the bases of the Global Agro-Ecological Zones methodology 
from FAO (Fischer et al., 2012), in order to estimate the agro- 
climatic attainable yield. 


© River basin 
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æ 5000 — 7000 
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Fig. 4 Net water provision for crops irrigation (m? ha~') in mainland Spain according to the homogenization process conducted 


from data in Hydrological Plans of Spain. 
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Thus, the yield constraints considered in our work were as 
follows: water deficit (from water requirements and irrigation 
water provision previously assessed) and soil salinity (in accor- 
dance with the estimated ECe and the salinity tolerance of 
giant reed). 

Regarding soil salinity, data of ECe and giant reed yields 
were gathered from the literature (Williams & Biswas, 2009 [In 
Singh, 2013]; Zema et al., 2012; Cosentino et al., 2013); after- 
wards, the respective decline in yield was calculated (Table 2). 

The yield constraint related to water deficit was assessed 
from an experiment conducted during four consecutive years 
in Sicily (Cosentino et al., 2014). Yield losses at the three differ- 
ent levels of water deficit tested — in terms of maximum evapo- 
transpiration restitution — were compiled (Table 3). The 
irrigation amount was inferred from the water-use efficiency 
and rainfall values. 

Hence, a function of yield loss was built for each parameter 
from literature data using 28.9 t dm ha~' (Cosentino et al., 
2014) as a reference value for giant reed potential yield, includ- 
ing stems and leaves. In this way, the agro-climatic attainable 
yield was estimated for the limiting conditions of the identified 
saline areas. 

Nevertheless, the above-mentioned reference value consid- 
ered as the optimal yield may still be rather optimistic for regu- 
lar conditions of water availability and soil quality in 
Mediterranean areas. Therefore, a sensitivity analysis was also 
conducted by reducing this reference value by 20% and 50% 
(23.1 and 14.4 t dm ha“, respectively), to which the yield loss 
functions were then applied. Thus, the lignocellulosic biomass 
potential production from giant reed in standard conditions 
under water and salinity stress was also estimated. 


Results 


Identification of saline and saline-prone areas 


Of the 2 198 264 ha of the permanently irrigated agricul- 
tural land in Spain mapped by CLC2006, the area that 
was located within the saline soils area identified by Tóth 
et al. (2008) amounted to 644 581 ha. The irrigated 
agricultural land located in saline soils, broken down by 
origin of the irrigation water, was 25 680 ha for ground- 


Table 2 Yield losses of giant reed vs. electrical conductivity 
(ECe) calculated from mean values reported by the respective 
authors. Significance of the relationship between ECe and yield 
loss: P < 0.001 


Source ECe (dS m“!) Yield loss (%) 
Cosentino et al. (2013) 2,2 0.0 
6.3 15.7 
9.1 44.5 
Zema et al. (2012) 0.38 0.0 
0.6 2.4 
Williams & Biswas 8.3 25.0 
(2009) [In Singh, 2013] 12 50.0 
25 100.0 


Table 3 Data of yield losses of giant reed vs. water deficit, 
calculated from mean values reported by Cosentino et al., 2014. 
Significance of the relationship between water deficit and yield 
loss: P < 0.001 


Year of experiment Water deficit (%) Yield loss (%) 


1998/1999 0 0.0 
50 13.7 
70 29.8 

1999/2000 0 0.0 
50 12.2 
74 34.6 

2000/2001 0 0 
50 21.9 
76 40 


water, 141 474 ha for surface water and 477 427 for 
undifferentiated origin-type water, according to the 
information compiled from the Hydrological Plans. 

Regarding the interpolation of ECw point data, the 
Inverse Distance Weighting (IDW) interpolation method 
was the one that best performed among the methods 
tested for both groundwater and surface water control 
networks, according to the root-mean-square error 
(RMSE). The validation of the interpolation from the 
integrated database (with the groundwater and surface 
networks combined) showed a RMSE of 8.744 dS m'! 
and a mean relative error of 30%. IDW data interpola- 
tion showed that ECw ranged from 0.192 to 
15.502 dS m™!. According to that, ECe values varied 
from 0.288 to 23.253 dS m'. 

As it can be seen in Fig. 5, soil salinity in most irri- 
gated agricultural lands of mainland Spain could be 
classified as low (ECe < 4 dS m“') or moderately low 
(4<ECe<6dSm_/!). The Segura river basin (XIII in 
Fig. 5) showed the highest proportion of land with esti- 
mated values of ECe over 6 dS m™'. River basins V 
(Eastern Cantabria) and VII (Galician coast) did not 
include irrigated land according to CLC2006. 

Thus, according to the analysis conducted in this 
work there are 65 588-ha saline-prone agricultural lands 
in mainland Spain with 6 < ECe < 9 dS m_ |, of which 
34 412 ha are located on the saline soils identified by 
Toth et al. (2008). This latter saline land area is mainly 
concentrated in the Segura River Basin (Fig. 6) and was 
taken in our study as the saline agricultural land area 
suitable for growing giant reed without food market 
distortion. 


Estimation of crop requirements and the resource 
availability 


The annual growth cycle of giant reed in the Spanish 
conditions was found to extend from early March to late 
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Fig. 5 Distribution of salinity classes of the irrigated agricultural land within each River Basin according to the electrical 
conductivity of the soil extract (ECe). Key to River Basins as in Fig. 2. 
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Fig. 6 Saline-irrigated agricultural land suitable for growing giant reed. 


November, depending on the region. Kc for the growth 
cycle was estimated per set of 10 days, but for better 
interpretation, it is displayed in Fig. 7 per fortnight. 

The irrigation water requirements for optimal yield of 
giant reed in Spain were estimated at the range 2500- 
11 700 m? ha™! yr! taking into account the values of 
Kc, PET and Pef during the annual growth cycle 
(Fig. 8). Specifically, in the saline-irrigated agricultural 
lands identified as suitable for growing giant reed 
(6<ECw<9dSm/'!)_ the irrigation requirements 
varied between 7918 and 10 710 m? ha™! yr. 


Water provision for crop irrigation established by the 
Hydrological Plans in the potential giant reed growing 
area varied from 1226 to 5585 m? ha! yr!, according 
to our approach. 


Biomass potential production 


As the identified saline area has limiting conditions for 
plant growth (salinity and water deficit), the effect of 
yield constraints was estimated from yield reduction 
functions built for giant reed. 
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Fig. 7 Map of the evolution of the crop coefficient (Kc) of giant reed during its growth cycle in mainland Spain as estimated per fort- 


night (starting day—ending day/month). 


For soil salinity constraint, the best fit (R? = 0.970) 
from compiled data (Table 2) was found for a linear 
function as shown in Fig. 9. It was expressed as (7): 


Yi sai(%) = 0.042 x ECe (dS m™') — 0.0386 (7) 


where Yjsa1 is the yield loss (%) in relation to the 
potential crop yield due to salinity conditions and ECe 
is the electrical conductivity of the soil extract in 
dS m™. 

Regarding the water-deficit constraint, a polynomial 
equation (R? = 0.966) was obtained (Fig. 10) from data 
shown in Table 3. The equation is given in (8): 


Yiwat(%) = 0.71 x WD?(%) — 0.0437 x WD(%) + 0.0012 
(8) 


where Yiwat is the yield loss (%) due to water deficit 
and WD is the water deficit measured in percentage of 
the estimated irrigation water requirement. 

The yield reduction of giant reed in the identified sal- 
ine area ranged from 21.5% to 33.2% and from 1% to 
28%, for the salinity and water-deficit limiting condi- 
tions, respectively. Therefore, the agronomically attain- 
able yield considering these limiting growing conditions 
resulted in the range 12.7—22.2 t dm ha! yr! (Fig. 11). 
Hence, the potential production of lignocellulosic 


biomass from giant reed under the limiting growth con- 
ditions of the identified saline areas was estimated at 
597 338 t dm yr’. 

The sensitivity analysis resulted from the reduction of 
the assumed optimal yield in 20% and 50% (23.1 and 
14.4 t dm ha™', respectively) showed agronomically 
attainable yields with minimum values between 
10.2 t dm ha™ (20% reduction) and 6.4 t dm ha! (50% 
reduction) and maximum values of 17.8 and 
13.1 t dm ha™!, respectively. The consequently biomass 
potential production would decrease to 477 770 t 
dm yr! and 298 564 t dm yr, respectively. 


Discussion 


Marginal lands have been defined by Soldatos (2015) as 
lands not worth cultivating with agro-food crops 
because of a number of bio-physical or economic con- 
straints such as low soil quality, water and salinity 
stress, extreme climate conditions, soil and terrain hand- 
icaps, long distance to the market or state intervention, 
among others. Hence, saline lands and saline-water irri- 
gated lands in which traditional food crops would not 
be feasible may be considered as marginal land, where 
no competition or displacement of food production 
would be caused. To this end, this study addressed the 
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Fig. 8 Continuous data (raster) of irrigation water requirements for optimal giant reed yield in mainland Spain. Nonsuitable areas 
are referred to regions where the growth cycle would not be completed. 
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Fig. 9 Yield reduction function (%) based on the electrical 
conductivity of the soil extract (dS m7’). 
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Fig. 10 Yield reduction function (%) based on water deficit 
(%) according to irrigation water requirement of giant reed. 


identification of that type of lands and the estimation of 
the potential production of lignocellulosic biomass from 
giant reed, a saline-tolerant species. 

The identification of marginal land to be devoted to 
the production of lignocellulosic biomass from giant 
reed pursues a number of positive effects in terms of 


soil restoration and phytoremediation: unlike annual 
crops, perennials provide larger quantities of residual 
organic matter that improve the buffer capacity of the 
soil against rapid changes in salinity (Darwish et al., 
2005; Acosta et al., 2011). Additionally, salinization of 
most soils in the south-east of Spain has been often 
induced by the use of poor-quality irrigation water. The 
introduction of drip irrigation systems would counteract 
this effect (Acosta et al., 2011). Drip irrigation can be 
used for growing giant reed instead of sprinkler or 
flood irrigation, systems that are often used for growing 
field crops in dry environments. Furthermore, cultiva- 
tion of salt-tolerant and high water-use-efficiency crops 
would contribute to decrease the risk of soil salinization. 
Soil erosion and polluted runoff are also reduced by the 
substitution of conventional crops grown on marginal 
lands for bioenergy with low-input perennial grasses 
(Campbell et al., 2008). Consequently, the work here 
presented could help policy makers to know the poten- 
tial size of the bioenergy resources while avoiding the 
abandonment of agricultural land and its progressive 
degradation. 

Toth et al. (2008) updated the map of salt-affected 
soils of Europe from two databases: i) the European Soil 
Database at 1 : 1 000 000 scale (ESDB), which character- 
izes distinct soil types grouped in soils associations, and 
ii) the map of salt-affected soils in Europe (Szabolcs, 
1974) which was an extension of the FAO/UNESCO 
Soil map of the World specifying salt limitations for 
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Fig. 11 Agronomically attainable yield (t dm ha~! yr~') of giant reed in the identified saline agricultural areas. 


agricultural activities. As it was not technically feasible 
to spatially delineate every characterized soil or soil 
association — Soil Typological Units (STU), they were 
grouped into Soil Map Units (SMU) which could be dig- 
itized, managed and depicted by a GIS. Therefore, the 
methodology used for mapping saline soils only 
allowed classifying SMUs according to the share of area 
covered by saline soils (more or <50% — see Fig. 1). 
Based on the valuable work by Tóth et al. (2008), the 
procedure followed in our study for downscale saline 
lands was focused on the estimation and identification 
of the saline or brackish water. Besides, this approach 
also allowed the identification of areas with saline- 
prone soils. 

Soil salinity could evolve over time either positively 
or negatively, as a function of natural processes (pri- 
mary salinization) and human interventions (secondary 
salinization). Taking into account the date of the publi- 
cation considered in this work as a reference for saline 
soils identification (the year 2008), the use of comple- 
ment databases to update soil base information would 
be adequate. To this end, the automatic and manual net- 
works of ECw measurements, that compile data in the 
rivers and water bodies of Spain yearly, were used. In 
any case, further saline soil information would be 
needed to improve the current analyses. 

Mapping soil salinity encounters difficulties due to 
large spatial and temporal variability of the saline 


condition. Many studies — aiming at the elaboration of 
saline soils maps — are based on survey instruments 
such as electromagnetic induction metres (Lesch et al., 
2005), which makes the mapping of large areas time- 
consuming and expensive. To lower survey costs, field 
spectroscopy and remote sensing have been widely 
used, but satellite imagery also presents limitations 
derived from spectral and spatial resolution (Fernan- 
dez-Buces et al., 2006). Thus, studies using remote sens- 
ing for this purpose have been applied to limited study 
areas (Kauter et al., 2003; Fernandez-Buces et al., 2006; 
Madyaka, 2008; Melendez-Pastor et al., 2012). On the 
contrary, our approach addresses a first approximation 
of mapping soil salinity areas at a national scale; besides 
the method here developed could be applied to other 
regions or study areas where similar geo-spatial infor- 
mation were available. 

In the study by Melendez-Pastor et al. (2012), different 
digital image processing techniques were applied on 
multispectral images in order to analyse soil salinization 
processes in a coastal zone of south-east Spain, which is 
mainly devoted to irrigated farming. The study area 
was found to present soil electrical conductivities (ECe) 
from 1.38 to 4.71 dS m™!. Even though the former work 
is not comparable to the present study in terms of work- 
ing scale and techniques, the ECe values estimated in 
our work for the same area ranged between 3.1 and 
6.9 dS m™!. Another study (Herrero & Pérez-Coveta, 
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2005) analysed the trends on salinization in the Flumen 
irrigation district of Aragon (Northeast Spain) by direct 
measurements of soil salinity from three surveys (1975, 
1985/1986 and 1999). The median ECe surveyed in 1999 
for the considered saline soils was 1.94 dS m™!, a simi- 
lar value to the ECe values estimated in our study for 
that study area, which ranged between 1.3 and 
1.9 dS m |. Therefore it can be inferred an acceptable 
correlation between the results of the above-mentioned 
studies and our values of estimated ECe. Nevertheless, 
the saline land areas identified in the current paper 
should be validated by measuring the ECe of soil sam- 
ples in further studies. 

Regarding the potential area for giant reed in Spain 
— identified according to crop requirements — the 
results obtained in our work are consistent with the 
distribution map of possible energy crops (including 
giant reed) in Europe reported Zegada-Lizarazu et al. 
(2010). Both studies show that the distribution area of 
giant reed covers most mainland of Spain, except for 
the highlands of NE and NW Spain where the period 
of time that would be suitable for growing giant reed 
is shorter than the duration of the growth cycle of 
this crop. 

The Kc values estimated in this study are slightly 
higher than those adopted by Cosentino ef al. (2014), 
who assumed the values they used for Miscanthus x 
giganteus (Cosentino et al., 2007) to calculate the daily 
ETc and the irrigation needs. The irrigation crop 
requirements found in our work for Spain (2500— 
11 700 m? ha“! yr') were in line with the values 
reported for this crop in other studies, that is from 7360 
to 12 350 m? ha! in the Mediterranean conditions of 
Italy (Borin et al., 2013); 10 230 m° ha! under optimum 
water supply conditions in southern Italy (Zegada- 
Lizarazu et al., 2010) and 7000 m? ha! for semi-arid 
conditions in Greece (Dalianis et al., 1995). 

In our work, the potential saline area for growing 
giant reed was restricted to the ECe range of 6- 
9 dS m | (4-6 dS m | ECw) in order to conduct a con- 
servative assessment. As stated in other studies, giant 
reed may tolerate higher values of electrical conductiv- 
ity like 12 dS m™! ECe (Williams & Biswas, 2009) or 
16 dS m` ECw (Sánchez et al., 2015). If such values 
were assumed for our study, the potential area for 
growing giant reed in saline lands of Spain would be 
up to 84 885 ha and 87 782 ha, respectively. In any case, 
irrigation water with salinity of ECw > 2.4 dS m ! is 
highly saline and considered not suitable for irrigation 
under ordinary conditions for traditional field crops 
(e.g. cereals and citrus) (Williams & Biswas, 2010). Thus, 
there would not be a considerable displacement of food 
production by the devotion of the identified saline area 
to giant reed cultivation. 


In this regard, research has been recently conducted 
with the aim of identifying giant reed clones tolerant to 
water and salinity stress. In the study by Sanchez et al. 
(2015), the clones ‘Martinensis’ and ‘Piccoplant’ per- 
formed best in both limiting conditions. According to 
that study, giant reed seems to be more tolerant to salin- 
ity than to water deficit. Nevertheless, further develop- 
ment of highly salinity- and drought-tolerant clones of 
giant reed would be recommendable. 

Limiting growing conditions that combine the factors: 
water deficit and saline conditions were assessed as the 
sum of both yield reduction factors for this study. In the 
study by Sanchez et al. (2015), the yield reduction 
observed in the treatment that combined water stress 
(25% field capacity) and saline conditions (16 dS m~! 
ECw) was 66.3%, which represents a lower figure than 
for the sum of the separated limiting factors — 61.5% for 
water stress and 45.8% for salinity. In our work, the area 
identified with the harshest conditions — in terms of 
water availability and soil salinity — combined 64% 
water deficit and 7.9 dS m~! ECe, which led to 26.4% 
and 29.5% yield reduction, respectively. From that, the 
agronomically attainable yield was estimated at 
12.7 t dm ha! yr | for those conditions, which repre- 
sented 55.9% combined yield reduction, leading to a 
more conservative assessment than the results obtained 
by Sanchez et al. (2015). Studies on the separate effect of 
water stress and saline conditions on biomass yield can 
be found in the literature (Williams & Biswas, 2010; 
Cosentino et al., 2013, 2014). Nonetheless, to the best of 
our knowledge, there is no further research on the com- 
bining effect of both limiting factors when acting 
together, apart from the study by Sánchez et al. (2015). 
It is worth noting that the effect of stress factors on bio- 
mass yield in that experiment was studied for a period 
of only 2 months, not achieving the completion of the 
growth cycle of giant reed. In addition, the experiment 
was conducted with pot-grown plants. Hence, the effect 
of concurrent conditions of salinity and drought on 
giant reed yield needs further investigation. 

Taking the factor of soil salinity separately, Sanchez 
et al. (2015) found statistical differences between the 
ECw treatment of 16 dS m™! and the control treatment 
for total plant dry weight but not between 4 dS m ! 
and the control. According to the yield reduction func- 
tion estimated in our work, those conditions of salinity 
would lead to 96.6% and 21.3% yield reduction, respec- 
tively. 

The methodology developed for the estimation of 
actual yields in the present study, based on yield limit- 
ing factors, was somewhat similar to the analysis con- 
ducted by Wicke et al. (2011). Their assessment of the 
technical and economic potential of forestry plantations 
on salt-affected soils, at a global scale, was based on 
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limiting factors applied on constraint-free yield of tree 
species, while saline soils were determined from the 
Harmonized World Soil Database (HWSD) (FAO et al., 
2008). Several factors hamper the comparison between 
both studies, namely the spatial resolution and scale of 
the analysis, irrigation assumptions, the energy crops 
considered and the land-use databases used, among 
others. Thus, the actual yields estimated by Wicke et al. 
(2011) ranged between 0 and 27 t dm ha! yr! with 
3.1 t dm ha! yr! mean yield, whereas in our estima- 
tions the agronomically attainable yield of giant reed 
resulted in the range 12.7-22.2 t dm ha! yr (17.4 t 
dm ha! yr! weighed mean). Besides, Wicke et al. 
(2011) considered the whole range of values for soil 
salinity (i.e. ECe > 16.0 dS m™! but associated to low 
crop suitability and low yield) in their analysis. In con- 
trast, we removed from our analysis the agricultural 
land with ECe > 9 dS m1, even though giant reed 
would be able to survive in such conditions yielding 
very low amount of biomass. 

The potential biomass production from giant reed in 
the saline lands identified in this work resulted in 
597 338 t dm yr, an amount equivalent to 10.5-10° GJ 
yr! primary energy, assuming 17.58 MJ ke! low 
heating value (LHV) of (https://www.ecn.nl/phyllis2/ 
). Taking into account the yields assessed in the sensi- 
tivity analysis of the reference value (23.1 and 14.4 t 
dm ha™?), the potential primary energy production 
would result in 8.4-10° and 5.2-10° GJ yr, respec- 
tively. In this regard, it should be mentioned that the 
calorific value of giant reed biomass obtained from 
marginal soils might be different from the biomass 
produced in standard conditions. Several studies have 
been conducted focusing on the effect of water and 
salinity stress on yields, morpho-biometric and physio- 
logical parameters of energy crops (Sidella et al., 2015; 
Zanetti etal., 2015). However, to the best of our 
knowledge, the response of plants to salinity and 
drought in terms of the biomass energy value has not 
been deeply analysed. In any case, other factors 
regarding biomass harvesting and logistics (moisture 
content and biomass fractionation, among others) 
might have more relevant effects on the biomass 
energy content than the marginality of the production 
area. In fact, experiments of giant reed grown in 
water-deficit conditions in Madrid (Spain) carried out 
by the authors (publication pending) showed no statis- 
tical differences between water regimes for LHV. The 
weighed LHV mean was found 17.57 + 0.42 MJ kg! 
(2.4% variation coefficient), taking into account the bio- 
mass partitioning into stems (82.4 + 7.8%) and leaves 
(17.6 + 7.8%). This value is almost the same as the 
LVH assumed in our work. However, further research 
on this issue would be needed. 


The methodology developed in our approach can be 
applied to other study areas provided that geo-refer- 
enced data are available. Both methodology and results 
of the case study of Spain addressed the use of saline or 
saline-prone agricultural areas for biofuels avoiding 
competition with agro-food land uses while a crop cited 
as ‘non-food cellulosic material’ in the European Union 
Directive 2015/1513 is promoted. 


Acknowledgements 


The present work was supported by the European project 
OPTIMA (Optimization of Perennial Grasses for Biomass 
Production) http://www.optimafp7.eu/, 7th Framework Pro- 
gramme of the European Community for research, technologi- 
cal development and demonstration activities, Grant 
Agreement no. 289642. 


References 


Acosta JA, Faz A, Jansen B, Kalbitz K, Martinez-Martinez S (2011) Assessment of 
salinity status in intensively cultivated soils under semiarid climate, Murcia, SE 
Spain. Journal of Arid Environments, 75, 1056-1066. 

Allen RG, Pereira LS, Smith M, Raes D, Wright JL (2005) Dual crop coefficient 
method for estimating evaporation from soil and application extensions. Journal 
of Irrigation and Drainage Engineering, 131, 2-13. 

Araus JL, Slafer GA, Reynolds MP, Royo C (2002) Plant breeding and drought in C3 
cereals: what should we breed for? Annals of Botany, 89, 925-940. 

Ayers R, Westcot D (1985) Water quality for agriculture. In FAO Irrigation and Drai- 
nage Paper 29. Rev. 1. (ed. Food and Agriculture Organization of the United 
Nations). Food and Agriculture Organization of the United Nations, Viale delle 
Terme di Caracalla, Rome, Italy. Available at: http://www.fao.org/DOCReP/ 
003/T0234e/T0234e00.htm (accessed 15 July 2015) 

Borin M, Barbera AC, Milani M, Molari G, Zimbone SM, Toscano A (2013) Biomass 
production and N balance of giant reed (Arundo donax L.) under high water and 
N input in Mediterranean environments. European Journal of Agronomy, 51, 117— 
119. 

Bouaziz M, Matschullat J, Gloaguen R (2011) Improved remote sensing detection of 
soil salinity from a semi-arid climate in Northeast Brazil. Comptes Rendus - Geo- 
science, 343, 795-803. 

Calheiros CSC, Quitério PVB, Silva G, Crispim LFC, Brix H, Moura SC, Castro PML 
(2012) Use of constructed wetland systems with Arundo and Sarcocornia for pol- 
ishing high salinity tannery wastewater. Journal of Environmental Management, 95, 
66-71. 

Campbell JE, Lobell DB, Genova RC, Field CB (2008) The global potential of bioen- 
ergy on abandoned agriculture lands. Environmental Science and Technology, 42, 
5791-5794. 

Confederación Hidrográfica del Duero (2012) Project for the Hydrological Plan. Annex 
no. 5. Chapter no. 4. Water demands. 188 pp. 

Cosentino SL, Patanè C, Sanzone E, Copani V, Foti S (2007) Effects of soil water con- 
tent and nitrogen supply on the productivity of Miscanthus x giganteus Greef et 
Deu. in a Mediterranean environment. Industrial Crops and Products, 25, 75-88. 

Cosentino SL, Copani V, Testa G, Sidella S, Patanè G, Scordia D (2013) Giant Reed 
screening to salinity levels. In: 21st European Biomass Conference and Exhibition, 
Vol. 2, (ed. ETA-Florence Renewable Energies), pp. 147-150. ISBN 978-88-89407- 
53-0. 3-7 June 2013, Copenhagen, Denmark. 

Cosentino SL, Scordia D, Sanzone E, Testa G, Copani V (2014) Response of giant 
reed (Arundo donax L.) to nitrogen fertilization and soil water availability in semi- 
arid Mediterranean environment. European Journal of Agronomy, 60, 22-32. 

Curt MD, Barreiro I, Sanz M, Fernandez J (2014) Response of 1-node-cuttings of 
Arundo donax clones to salinity stress conditions. In: Abstract book Plant Biology 
Europe FESPB/EPSO 2014 Congress (ed. European Plant Science Organisation), pp. 
31. 22-26 June 2014, Dublin, Ireland. 

Dale VH, Kline KL, Wiens J, Fargione J (2010) Biofuels: implications for land use and 
biodiversity. Biofuels and Sustainability Reports. Ecological Society of America 
(ESA), Washington, DC. 


© 2015 The Authors. Global Change Biology Bioenergy Published by John Wiley & Sons Ltd., doi: 10.1111/gcbb.12329 


14 J. SANCHEZ et al. 


Dalianis C, Sooter C, Christou M (1995) Growth, biomass productivity of Arundo 
donax and Miscanthus sinensis “giganteus.” In: Biomass for Energy, Environment, 
Agriculture and Industry: Proceedings of the 8th EU Biomass Conference (ed. Al. C et), 
pp. 575-582. Pergamon Press, Vienna, Austria. 

Darwish T, Atallah T, Moujabber M, Khatib N (2005) Salinity evolution and crop 
response to secondary soil salinity in two agro-climatic zones in Lebanon. Agri- 
cultural Water Management, 78, 152-164. 

Dastane N (1978) Effective Rainfall in Irrigated Agriculture (ed. Food and Agriculture 
Organizario of Unated Nations). FAO Irrigation and Drainage paper 25. 

Di Nasso NN, Roncucci N, Triana F, Tozzini C, Bonari E (2011) Productivity of giant 
reed (Arundo donax L.) and miscanthus (Miscanthus x giganteus greef et deuter) as 
energy crops: growth analysis. Italian Journal of Agronomy, 6, 141-147. 

Douaoui AEK, Nicolas H, Walter C (2006) Detecting salinity hazards within a semi- 
arid context by means of combining soil and remote-sensing data. Geoderma, 134, 
217-230. 

Emberger L (1955) Une classification biogéographique des climats. Recueil des Tra- 
vaux du Laboratoire de Botanique et Géologie. Faculté des Sciences de Montpel- 
lier. Sér. Bot., 7, 3-43. 

European Commission and the European Soil Bureau Network (2004) The European 
Soil Database distribution version 2.0 (ESDB, 2004). 

FAO,IIASA,ISRIC,ISS-CSA,JRC (2008) Harmonized World Soil Database (version 1.0). 
FAO, Rome and IIASA, Laxenburg, Austria. 

Fernandez-Buces N, Siebe C, Cram S, Palacio JL (2006) Mapping soil salinity using a 
combined spectral response index for bare soil and vegetation: a case study in 
the former lake Texcoco, Mexico. Journal of Arid Environments, 65, 644-667. 

Fischer G, Nachtergaele F, Prieler S et al. (2012) Global Agro-Ecological Zones (GAEZ 
v3.0) — Model Documentation. 179 pp. 

Gommes R (1983) Pocket Computers in Agrometeorology. Food and Agriculture Organi- 
zation (FAO). Paper 45. Rome, Italy. 

Hargreaves GH (1994) Defining and Using Reference Evapotranspiration. Journal of 
Irrigation and Drainage Engineering, 120, 1132-1139. 

Hasanuzzaman M, Nahar K, Alam MM et al. (2014) Potential use of halophytes to 
remediate saline soils. BioMed Research International, 2014, 1-12. 

Herrero J, Pérez-Coveta O (2005) Soil salinity changes over 24 years in a Mediter- 
ranean irrigated district. Geoderma, 125, 287-308. 

Kang S, Post W, Wang D, Nichols J, Bandaru V, West T (2013) Hierarchical marginal 
land assessment for land use planning. Land Use Policy, 30, 106-113. 

Kauter D, Lewandowski I, Claupein W (2003) 03/02394 Quantity and quality of har- 
vestable biomass from Populus short rotation coppice for solid fuel use - a review 
of the physiological basis and management influences Kauter, D. et al. Biomass 
and Bioenergy, 2003, 24, (6), 411-427. Fuel and Energy Abstracts, 44, 389. 

Lesch SM, Corwin DL, Robinson DA (2005) Apparent soil electrical conductivity 
mapping as an agricultural management tool in arid zone soils. Computers and 
Electronics in Agriculture, 46, 351-378. 

Lewis S, Kelly M (2014) Mapping the potential for biofuel production on marginal 
lands: differences in definitions, data and models across scales. ISPRS Interna- 
tional Journal of Geo-Information, 3, 430-459. 

Madyaka M (2008) Spatial Modelling and Prediction of Soil Salinization Using SaltMod 
in a GIS Environment. International Institute for Geo-Information Science and 
Earth Observation, Enschede, the Netherlands. 

McMaster GS, Wilhelm WW (1997) Growing degree-days: one equation, two inter- 
pretations. Agricultural and Forest Meteorology, 87, 291-300. 

Melendez-Pastor I, Hernandez EI, Navarro-Pedreno J, Gómez I (2012) Mapping soil 
salinization of agricultural coastal areas in southeast Spain. In: Remote Sensing 
Applications (ed. B. Escalante), pp. 117-140. InTech, Available from: http:// 
www.intechopen.com/books/remote-sensing-applications /mapping-soil-saliniza 
tion-of-agricultural-coastal-areas-in-southeast-spain. 

Milbrandt A, Overend R (2009) Assessment of Biomass Resources From Marginal Lands 
in APEC Economies. National Renewable Energy Laboratory (NREL), Golden, Col- 
orado. 

National Meteorological Agency of Spain (AEMET) (2010a) Raster layers (1 km pixel 
size) of mean monthly temperatures from 1971-2000 period in ESRI™ format. 

National Meteorological Agency of Spain (AEMET) (2010b) Raster layers (1 km pixel 
size) of mean minimum temperatures from 1971-2000 period in ESRI™ format. 

National Meteorological Agency of Spain (AEMET) (2010c) Raster layers (1 km pixel 
size) of monthly precipitation from 1971-2000 period in ESRI™ format. 

Nijsen M, Smeets E, Stehfest E, van Vuuren DP (2012) An evaluation of the global 
potential of bioenergy production on degraded lands. GCB Bioenergy, 4, 130-147. 

Ochieng GM, Ojo OI, Otieno FA, Mwaka B (2013) Use of remote sensing and geo- 
graphical information system (GIS) for salinity assessment of Vaal-Harts irriga- 
tion scheme, South Africa. Environmental Systems Research, 2, 4. 


Rhoades J, Kandiah A, Mashali A (1992) Chapter 4 - Water quality assessment. In: 
The Use of Saline Waters for Crop Production (ed. Food and Agriculture Organiza- 
tion of the United Nations), pp. 1-133. FAO Irrigation and Drainage paper 48, 
Rome, Italy. 

Sanchez E, Scordia D, Lino G, Arias C, Cosentino SL, Nogués S (2015) Salinity and 
Water Stress Effects on Biomass Production in Different Arundo donax L. Clones. 
BioEnergy Research, 8, 1461-1479. 

Shortall OK (2013) “Marginal land” for energy crops: exploring definitions and 
embedded assumptions. Energy Policy, 62, 19-27. 

Sidella S, Testa G, Scordia D, Cosentino SL (2015) Screening of giant reed clones 
under increasing salinity levels. In: Perennial Biomass Crops for a Resource-Con- 
strained World, p. 68 Book of Abstracts. Poster — P4.05. 7-10 Sept. 2015. 

Singh BP (2013) Biofuel Crops: Production, Physiology and Genetics, Singh edn. B.P edn. 
CABI, Fort Valley State University, Fort Valley, Georgia. 

Soldatos P (2015) Economic aspects of growing perennial grasses for fibre and bioen- 
ergy in South Europe. In: Perennial Biomass Crops for a Resource-Constrained World, 
p. Oral Presentation — 02.06. 7-10 Sept. 2015, Stuttgart, Berlin. 

Szabolcs I (1974) Salt Affected Soils in Europe (ed. Nijhoff M) Research Institute for 
Soil Science and Agricultural Chemistry of the Hungarian Academy of Science, 
Budapest, Hungary. pp. 63. 

Tenerelli P, Carver S (2012) Multi-criteria, multi-objective and uncertainty analysis 
for agro-energy spatial modelling. Applied Geography, 32, 724-736. 

Toth G, Adhikari K, Varallyay G, Toth T, Bódis K, Stolbovoy V (2008) Updated map 
of salt affected soils in the European Union. In: Threats to Soil Quality in Europe 
(eds. Téth G, Montanarella L and Rusco E), EUR 23438 EN - 2008. pp. 65-77. 
European Commission. Joint Research Centre. Institute for Environment and Sus- 
tainability, Luxembourg. 

Trabucco A, Zomer RJ, Bossio DA, van Straaten O, Verchot LV (2008) Climate 
change mitigation through afforestation/reforestation: a global analysis of hydro- 
logic impacts with four case studies. Agriculture, Ecosystems and Environment, 126, 
81-97. 

Triana F, Nassi o Di Nasso N, Ragaglini G, Roncucci N, Bonari E (2014) Evapotran- 
spiration, crop coefficient and water use efficiency of giant reed (Arundo donax L.) 
and miscanthus (Miscanthus giganteus Greef et Deu.) in a Mediterranean environ- 
ment. GCB Bioenergy, 7, 811-819. 

Turley D, Taylor M, Laybourn R et al. , (2010). Assessment of the Availability of Mar- 
ginal and Idle Land for Bioenergy. Crop Production in England and Wales, Food 
and Environment Research Agency. Department for Environment Food and Rural 
Affairs (DEFRA). Research Project Final Report NF0444. London, UK. 

Wicke B, Smeets E, Dornburg V, Vashev B, Gaiser T, Turkenburg W, Faaij A (2011) 
The global technical and economic potential of bioenergy from salt-affected soils. 
Energy & Environmental Science, 4, 2669. 

Williams CMJ, Biswas TK (2009) Salinity tolerance, nutrient needs, weed risk man- 
agement and carbon sequestration of Arundo donax for remediation of high saline 
wastewaters for biofuel or pulp/paper feedstock production. Final Report Project 
No. 56/2008 for the South Australia Murray-Darling Basin Natural Resources 
Management Board, SARDI, Adelaide, Australia. [In Singh, 2013] 

Williams CMJ, Biswas TK (2010) Commercial Potential of Giant Reed for Pulp, Paper 
and Biofuel Production. Rural Industries Research and Development Corporation 
(RIRDC). Australian Government. RIRDC Publication No. 10/215. RIRDC Project 
No. PRJ-000070. 152 pp. ISBN 978-1-74254-180-8. 

Williams CMJ, Biswas TK, Black ID et al. (2009) Use of poor quality water to pro- 
duce high biomass yields of giant reed (Arundo donax L.) on marginal lands for 
biofuel or pulp/paper. International Symposium on Underutilized Plants for 
Food Security, Nutrition, Income and Sustainable Development, 806, 595-602. 

Zanetti F, Zegada-Lizarazu W, Sartoni R, Monti A (2015) Response of switchgrass 
cytotypes to salinity and temperature. In: Perennial Biomass Crops for a 
Resource-Constrained World, p. 69 Book of Abstracts. Poster — P4.06. 7-10 Sept. 
2015, Stuttgart, Berlin. 

Zegada-Lizarazu W, Elbersen HW, Cosentino SL, Zatta A, Alexopoulou E, Monti A 
(2010) Agronomic aspects of future energy crops in Europe. Biofuels, Bioproducts 
and Biorefining, 4, 674-691. 

Zema DA, Bombino G, Andiloro S, Zimbone SM (2012) Irrigation of energy crops 
with urban wastewater: effects on biomass yields, soils and heating values. Agri- 
cultural Water Management, 115, 55-65. 

Zhuang D, Jiang D, Liu L, Huang Y (2011) Assessment of bioenergy potential on 
marginal land in China. Renewable and Sustainable Energy Reviews, 15, 1050-1056. 
Zomer RJ, Trabucco A, Bossio DA, Verchot LV (2008) Climate change mitigation: a 
spatial analysis of global land suitability for clean development mechanism 
afforestation and reforestation. Agriculture, Ecosystems and Environment, 126, 67-80. 


© 2015 The Authors. Global Change Biology Bioenergy Published by John Wiley & Sons Ltd., doi: 10.1111/gcbb.12329 


